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In this paper we address the following question: do the recent advances in the orbit deter- 
mination of the major natural satellites of Saturn obtained with the analysis of the first 
data sets from the Cassini mission allow to detect the general relativistic gravitoelectric 
orbital precessions of the mean longitudes of such moons? The answer is still negative. 
The present-day down-track accuracy is adequate for Mimas, Enceladus, Thetys, Dione, 
Rhea and Titan and inadequate for Hyperion, Iapctus and Phoebe. Instead, the size of 
the systematic errors induced by the mismodelling in the key parameters of the Satur- 
nian gravitational field like the even zonal harmonics Ji are larger than the relativistic 
down-track shifts by about one order of magnitude, mainly for the inner satellites like 
Mimas, Enceladus, Thetys, Dione, Rhea, Titan and Hyperion. Iapetus and Phoebe are 
not sensibly affected by such kind of perturbations. Moreover, the bias due to the un- 
certainty in Saturn's GM is larger than the relativistic down-track effects for all such 
moons up to two orders of magnitude (Phoebe). Thus, it would be impossible to sep- 
arately analyze the mean longitudes of each satellite. Proposed linear combinations of 
the satellites' mean longitudes would allow to cancel out the impact of the mismodelling 
in the low-degree even zonal harmonics and GM, but the combined down-track errors 
would be larger than the combined relativistic signatures by a factor 10 3 . 

Keywords: Cassini spacecraft; general relativistic orbit precessions; Saturnian system of 
natural satellites 



1. Introduction 

A satisfactorily empirical corroboration of a fundamental theory requires that as 
many independent experiments as possible are conducted by different scientists in 
different laboratories. Now, general relativity is difficult to test, especially in the 
weak-field and slow-motion approximation, valid, e.g., in our Solar System, both 
because the relativistic effects are very small and the competing classical signals 
are often quite larger. Until now, Solar System tests of general relativity accurate 
to better than 1% have been performed by many independent groups only in the 
gravitational field of the Sun by checking the effects induced by the gravitoelectric 
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Schwarzschild a part of the space-time metric on the propagation of electromagnetic 
waves and planetary orbital motions 4 . It is as if many independent experiments 
aimed to measure fundamental physical effects were conducted always in the same 
laboratory. Thus, it is worthwhile to try to use different laboratories, i.e. other grav- 
itational fields, to perform such tests, even if their outcomes should be less accurate 
than those conducted in the Sun's field. In principle, the best candidates other than 
the Sun are the giant planets of the Solar System like Jupiter and Saturn. To date, 
the only investigations of this kind are due to Hiscock and Lindblom 5 , who prelim- 
inarily analyzed the possibility of measuring the Einstein pericentre precessions 6 of 
some of the natural satellites of Jupiter and Saturn, and to Iorio and Lainey 7 who 
investigated the measurability of the Lense-Thirring precessions in the system of 
the galilean satellites of Jupiter with modern data sets b . 

The aim of this paper is to investigate if the recent improvements in the 
ephemerides of the major Saturnian satellites (see Table 1 for their orbital pa- 
rameters) obtained from the first Cassini data allows to detect at least the general 
relativistic gravitoelectric precessions of their orbits. 

Table 1. Orbital parameters of the major 
Saturnian satellites: a is the semimajor axis, 
e is the eccentricity and i is the inclina- 
tion, referred to the local Laplace planes 
(http://ssd.jpl.nasa.gov/7sat_elcm). 



Satellite 


a (km) 


e 


i (deg) 


Mimas 


185540 


0.0196 


1.572 


Enccladus 


238040 


0.0047 


0.009 


Thetys 


294670 


0.0001 


1.091 


Dionc 


377420 


0.0022 


0.028 


Rhea 


527070 


0.0010 


0.331 


Titan 


1221870 


0.0288 


0.280 


Hyperion 


1500880 


0.0274 


0.630 


Iapetus 


3560840 


0.0283 


7.489 


Phoebe 


12947780 


0.1635 


175.986 



2. The relativistic effects investigated 

The Einstein pericentre precession is undoubtedly the most famous relativistic or- 
bital effect, but it is not the only one. Indeed, also the mean anomaly M. undergoes 



a In regard to the much smaller gravitomagnetic Lense-Thirring precessions of the planetary orbits, 
they lie just at the edge of the present-day accuracy in planetary ephemerides and have recently 
been found in agreement with the latest measurements, although the errors are still large 2 . A 
6% test performed in the gravitational field of Mars with the MGS spacecraft has recently been 
reported in Ref. 3. 

b In their original papers Lense and Thirring 1 proposed to use some of the satellites of Jupiter and 
Saturn. 
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a secular advance which is even larger than that experienced by the pericentre. For 
small eccentricities it is 8 

9nGM m 
M GE — ===, (1) 



c 2 aVT 



where G is the Newtonian constant of gravitation, M is the mass of the central 
body which acts as source of the gravitational field, a, e and n = \J GM/a 3 are 
the semimajor axis, the eccentricity and the Keplerian mean motion, respectively, 
of the satellite's orbit. Eq. (1), which yields a rate of -130 arcseconds per century 
for Mercury, has been obtained by using the the standard isotropic radial coordi- 
nate r related to the Schwarzschild coordinate r by r = r(l + GM/2c 2 r) 2 . The 
validity of Eq. (1) has also been numerically checked by integrating over 200 years 
the Jet Propulsion Laboratory (JPL) equations of motion of all the planets of the 
Solar System, which are written in terms of just the standard isotropic radial co- 
ordinate, with and without the gravitoelectric 1/c 2 terms in the dynamical force 
models 9 in order to single out just the post-Newtonian gravitoelectric effects (E.M. 
Standish, private communication, 2004). The obtained precessions fully agree with 
those obtainable from Eq. (1). 

As a consequence, the mean longitude A = u> + £1 + A4, which is one of the non- 
singular orbital elements used for orbits with small eccentricities and inclinations 
like those of most of the Solar System major bodies as just the Saturnian satellites 
(see Table 1), experiences a secular advance 

Age ~ 5 — ■ (2) 

It is twice the pericentre advance. On the other hand, because of M. in the definition 
of A the systematic uncertainty in the Keplerian mean motion n must also be 
accounted for if a secular rate must be extracted from the analysis of such an 
orbital clement. 



3. The possibilities offered by the Saturnian system 

Is it possible to measure the gravitoelectric orbital precessions of the natural satel- 
lites of Saturn, in view of the recent refinements of their ephemerides obtained by 
analyzing the first data from the Cassini spacecraft? In this Section we will ad- 
dress this problem in detail. To this aim, we must confront the magnitude of the 
relativistic effects of interest with the major systematic errors induced by classical 
forces having the same signatures (Section 3.1), and with the currently available 
measurement accuracies (Section 3.2). 

3.1. The reduction of the impact of the even zonal harmonics 

In regard to the first issue, a major source of systematic bias is represented by the 
even zonal harmonic coefficients Je, I = 2, 4, 6, ... of the multipolar expansion of the 



c Here ui and Q are the argument of pericentre and the longitude of the ascending node, respectively. 
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Saturn's gravitational potential. Indeed, they induce secular precessions on A 

Aeven zonals = ^ ' ^.i^ii (3) 

where the coefficients \j are expressed in terms of the satellite's orbital elements 
a, e, i and the mass and the radius of the central body. Such classical advances must 
be accurately modelled because they are much larger than the relativistic ones. 
Even the latest determinations of the Saturnian gravity field from some Cassini 
data sets 10 show that the currently available models of the even zonal harmonics 
do not reach the required accuracy to allow for a measurement of the Einstein 
precessions. This fact is clearly shown by Table 2. Indeed, we can note that for 

Table 2. Mean longitudes of the major Saturnian satel- 
lites: general relativistic gravitoelectric secular preces- 
sions AqEi in " cy , mismodcllcd classical secular pre- 
cessions <5A j 2 , i5A j 4 , in " cy~ 1 , due to the uncertainties in 
the even zonal harmonics J2 , J4 , systematic errors Sugm, 
in " cy -1 , in the Keplerian mean motions n due to the 
uncertainty in Saturn's GM. The values 8(GM) = 1.2 
km 3 s~ 2 , 8J 2 = 0.4 X 1(T 6 , 6J4 = 3 X lCT 6 have been 
used. 



Satellite 


Age 


<5Aj 2 


SXj 4 


SncM 


Mimas 


-684.670 


5926.552 


5465.353 


793.461 


Enccladus 


-367.207 


2478.985 


1390.913 


546.018 


Thetys 


-215.374 


1173.870 


429.231 


396.440 


Dionc 


-116.003 


493.910 


110.229 


273.492 


Rhea 


-50.334 


153.451 


17.558 


165.721 


Titan 


-6.152 


8.101 


0.172 


46.950 


Hyperion 


-3.679 


3.942 


0.055 


34.487 


Iapctus 


-0.424 


0.186 


e>(io- 4 ) 


9.437 


Phoebe 


-0.016 


0.002 


e>(io- 7 ) 


1.361 



Mimas, Enceladus and Thetys the mismodcllcd precessions due to J2 and J4 are 
larger than the relativistic rates by about one order of magnitude. For Dione, Rhea, 
Titan and Hyperion the bias due to J4 is smaller than the relativistic signal. For 
Iapetus and Phoebe the even zonal harmonics of Saturn do not represent a problem. 
In addition to the bias due to the {Ji}, for the mean longitude there is also the 
systematic error 5n GM — 5 (GM) / V AGMa 3 induced on the Keplerian mean motion 
n by the uncertainty in the Saturn's GM which, as can be seen, is larger than the 
relativistic shifts for all the satellites. 

A way to overcome the problem of the impact of the even zonals is represented 
by the linear combination approach 8 ' 11 summarized in the following. Let us write 
down the expressions of the measured residuals d of the mean longitudes <5A moas of 



The Keplerian orbital elements are not directly observable, so here we are using the words 'mea- 
sured residuals' in an improper sense. We mean, instead, a set of solve-for parameters of a suitable 
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N chosen satellites in terms of the mismodelled Keplerian mean motions, of the 
classical precessions due to the mismodelled parts of the even zonal harmonics, and 
of the relativistic precessions, assumed as totally unmodelled features of motion 

2(N-2) 

S*£L = n. GM 5(GM) + ]T A ( f SJ t + Ag^, 3 = b 2..JV, (4) 

e>2 

where X.e are the coefficients of the classical precessions of degree I of the mean 
longitude and 

n gm = , (5) 
VAGMa 3 

If we look at Eq. (4) as a inhomogeneous algebraic system of N linear equations 
in the N unknowns {5(GM); 5Je; £} and solve it for £, we get a linear combination 
of the satellites' mean longitudes residuals which, by construction, is independent 
of 5(GM) and the first N — 2 even zonal harmonics, and is sensitive just to the 
relativistic signatures and to the perturbations of the remaining, uncancelled even 
zonal harmonics. The cancellation of the low-degree even zonals is important also 
because in this way one avoids any possible a priori 'imprint' of the relativistic effects 
themselves via such spherical harmonics. Indeed, the Saturnian gravity field models 
are obtained as least-square solutions using all the available data from the spacecraft 
encounters and the satellites' motions around Saturn, so that the relativistic effects 
themselves are in some way included in the solved- for parameters like { Jg\. 
A possible combination which cancels out 8J2, SJ4 and 8{GM) is 

£ _ ^mcas + fcl^^mcas + ^2^A^ as + fc 3 (5A^' eas ^ 
^GE + ^l^GE "b ^2A G g + k^X^ 

with 

' fci = -0.1296, 
k 2 = 0.0224, 

(7) 

k 3 = -0.0067, 

, ^ge + ^ge + *2<^. + fc 3 ATi E = 0.0237 "cy- 1 

As can be argued from Table 2, the impact of the remaining uncancelled even zonal 
harmonics Jq, Jg is negligible. 



least-square solution to be obtained by contrasting all the available observational data to a set of 
dynamical force models in which all the 1/c 2 terms have purposely been switched off. 
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3.2. The present-day orbital accuracy 

In regard to the feasibility of such a measurement, it is useful to consider the 
down-track shifts At <~ aAA and the present-day accuracy according to the latest 
Saturnian ephemerides SAT240 12 over 81 years. The results arc in Table 3. It can 

Table 3. Major Saturnian satellites: grav- 
itoelectric down-track shifts, in km, and 
present-day orbital accuracy, in km, over a 
time span of 81 years. See also on the Inter- 
net http:/ /ssd.jpl.nasa.gov/?sat_ephcm. 



Satellite 


Atge (km) 


<5(Ar) (km) 


Mimas 


-498 


50 


Enccladus 


-343 


20 


Thetys 


-249 


50 


Dionc 


-171 


20 


Rhea 


-104 


20 


Titan 


-29 


10 


Hyperion 


-21 


100 


Iapctus 


-6 


20 


Phoebe 


-0.8 


50 



be noted that the situation is presently not favorable just for the satellites used 
for the combination of Eq. (6). This fact reflects on Eq. (6) itself: indeed the total 
uncertainty calculated by summing in quadrature the errors 5 (At) of Table 3 with 
the coefficients (7) 



^/[S(AT) Ph ] 2 + [M(Ar) Ia ] 2 + M(Ar) H y] 2 + [M(AT) Ti ] 2 



(8) 



At-P£ + fa At£ e + k 2 Ar|£ + k 3 At£ e 

is <~ 10 3 times larger than the gravitoelectric shift over 81 years. 

About the use of the inner satellites, for which the relativistic shifts Atqe are 
larger than the errors 5 (At), Table 2 shows that it is not possible to use their mean 
longitudes without combining them because both the even zonal mismodelled shifts 
and the bias of n are quite larger. On the other hand, it turns out that also the linear 
combination approach does not yield good results. Indeed, the root-sum-square of 
the errors of a combination with Enceladus, Dione, Rhea and Titan is 10 3 times 
larger of the relativistic combined shifts. 

4. Discussion and conclusions 

In regard to the measurability of the general relativistic gravitoelectric orbit ad- 
vances in the Saturn's system of natural satellites, it turns out that the present-day 
improvements in their ephemerides by the Cassini spacecraft are not yet sufficient 
to detect such post-Newtonian effects. 

As a general rule, in regard to the single satellites, an about one order of magni- 
tude improvement in the knowledge of the parameters of the Saturnian gravity field 
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like the even zonal harmonics Ji and GM would be required to make the compet- 
ing classical bias smaller than the relativistic precessions. This would make possible 
to successfully analyze separately the mean longitudes of each of Mimas, Ence- 
ladus, Thetys, Dione, Rhea and Titan because the present-day down-track errors 
are already smaller than the relativistic shifts. The use of the mean longitudes of 
Hyperion, Iapetus and Phoebe, for which the Saturnian even zonals represent a rel- 
atively less important problem, is strictly related to a one-two orders of magnitude 
improvements in the down-track parts of their orbits. 

The linear combination approach yields, in principle, good results in reducing 
the impact of the mismodelling in the even zonals and GM. The main problem 
with such strategy relies in the down-track orbital accuracy. Indeed, the down-track 
errors weighted by their combination's coefficients must be summed in quadrature, 
while the weighted relativistic shifts are combined with their proper own signs, so 
that with the present-day precision the overall error amounts to 10 3 . 
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